In [24]:
%matplotlib inline
Vamos elaborar um diagrama TS com o auxílio do pacote gsw [https://pypi.python.org/pypi/gsw/3.0.3], que é uma alternativa em python para a toolbox gsw do MATLAB:
In [25]:
import gsw
Se você não conseguiu importar a biblioteca acima, precisa instalar o módulo gsw.
Em seguida, importamos a biblioteca numpy que nos permite usar algumas funções matemáticas no python:
In [26]:
import numpy as np
import matplotlib.pyplot as plt
In [28]:
sal = np.linspace(0, 42, 100)
temp = np.linspace(-2, 40, 100)
s, t = np.meshgrid(sal, temp)
# Abaixo usamos diretamente o resultado da biblioteca gsw:
# Thermodynamic Equation Of Seawater - 2010 (TEOS-10)
sigma = gsw.sigma0(s, t)
# Quantidade de linhas desejada
cnt = np.arange(-7, 35, 5)
fig, ax = plt.subplots(figsize=(5, 5))
ax.plot(sal, temp, 'ro')
# O comando abaixo faz curvas de nível com dados contour(X, Y, Z)
cs = ax.contour(s, t, sigma, colors='blue', levels=cnt)
# Aqui fazemos rótulos para as curvas de nível
ax.clabel(cs, fontsize=9, inline=1, fmt='%2i')
ax.set_xlabel('Salinity [g kg$^{-1}$]')
ax.set_ylabel('Temperature [$^{\circ}$C]')
#plt.plot(s,t,'ro')
Out[28]:
In [ ]: